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ABSTRACT 


The  reaction  32S(d,n)33Cl  has  been  investigated  by  the  neutron 
time-of-flight  method.  Precise  reaction  Q-values  and  the  correspon¬ 
ding  excitation  energies  for  the  low-lying  levels  of  33C1  have  been 
measured.  Angular  distributions  of  the  neutrons  leading  to  the 
ground,  0.812,  2.688  and  2.851  MeV  levels  have  been  studied  at 
Ej  =  4.70  and  5.50  MeV.  Angular  distributions  of  the  neutrons 
leading  to  the  1.993  and  2.351  MeV  states  have  been  measured  at 
=  4.54  and  4.70  MeV.  The  DWBA  theory  has  been  employed  to 
determine  the  orbital  angular  momentum  l  of  the  captured  proton. 

The  Hauser-Feshbach  calculations  have  been  performed  to  see  how 
important  is  the  contribution  from  the  compound  nucleus  formation. 
Reasonable  agreement  has  been  obtained  between  DWBA  theory  and 
experimental  results.  7-values  of  2,  0,  3  and  1  have  been  obtained 
for  the  ground,  0.812,  2.688  and  2.851  MeV  states,  respectively. 

The  absolute  spectroscopic  factors  have  been  extracted  and 
compared  to  shell  model  predictions.  Qualitative  agreement 
has  been  obtained. 
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CHAPTER  I 
INTRODUCTION 

Nuclear  reactions  constitute  one  of  the  most  powerful  tools 
of  nuclear  spectroscopy,  i.e.,  the  mapping  of  nuclear  energy  levels. 
The  level  schemes,  with  all  the  pertinent  information  about  each  state, 
of  course,  constitute  the  most  valuable  data  for  testing  theories  of 
nuclear  structure. 

The  present  work  reports  on  the  low-lying  states  of  33C1 
as  observed  in  the  32S(d,n)33Cl  reaction.  This  reaction  has  been 
studied  with  low  resolution  by  Middleton  ejt  a]_.  (Mi  53),  Macefield 
and  Towle  (Ma  60)  and  Mubarakmand  and  Macefield  (Mu  67).  In  a  (d,n) 
experiment  done  by  Grandy  et_  al_.  (Gr  68)  in  this  laboratory  on  a 
42Ca  target,  which  was  contaminated  with  32S,  a  number  of  transi¬ 
tions  to  states  in  33C1  have  been  identified.  The  levels  reported 
in  (d,n),  (3He,d)  and  (p,y)  experiments  on  32S  are  presented  in  fig. 
1-1.  It  is  noticed  that  the  levels  at  2.11  and  2.53  MeV  have  been 
seen  only  in  the  (d,n)  experiments,  while  the  levels  at  1.99  and  2.35 
MeV  have  been  observed  only  in  the  (3He,d)  and  (p,y)  work.  The  level 
at  2.69  MeV,  reported  in  (3He,d)  experiments,  also  has  been  seen  in 
Grandy  (d,n)  work  together  with  three  new  levels  at  2.63,  3.24  and 
3,27  MeV.  The  results  of  Grandy  et  al_.  are  somewhat  uncertain  since 
32S  was  just  a  contaminant. 


Figure  1-1.  Levels  of  the  mirror  nuclei  33C1  and  33S  as  observed 


in  different  reactions. 
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It  was  considered  that  a  study  of  the  32S(d,n)33Cl  reaction 
with  high  resolution  would  be  worthwhile  from  two  points  of  view:  to 
determine  precisely  the  excitation  energy  of  33C1  levels  excited  by 
the  (d,n)  reaction,  and  to  measure  the  angular  distribution  of  neutrons 
to  determine  the  angular  momentum,  parity  and  spectroscopic  factors  of 
the  33C1  levels. 

In  the  present  work,  time-of-f light  spectra  have  been  taken 
at  three  deuteron  energies  4.54,  4.70  and  5.50  MeV.  The  excitation 
energies  of  the  33C1  levels  have  been  measured.  Angular  distributions 
have  been  obtained  for  most  of  the  levels  and  compared  with  the  Dis¬ 
torted  Wave  Born  approximation  and  Hauser-Feshbach  theories.  Absolute 
spectroscopic  factors  have  been  extracted  and  compared  with  the  shell 
model  predictions. 


'  ‘  <  .  :  *r-  r  ■  n  o  ! 
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CHAPTER  II 

THE  EXPERIMENTAL  SETUP 
2.1  Beam  Transport  System 

To  date,  the  most  exact  work  in  neutron  spectroscopy  has 
undoubtedly  been  performed  with  pulsed  accelerators  and  electronic 
timing  devices  to  measure  the  time  of  flight  of  neutrons  over  a 
certain  path  length  from  the  target  to  the  detector.  Such  a  technique 
has  the  advantages  of  good  energy  resolution  and  applicability  to  a 
wide  range  of  neutron  energies. 

Since  in  the  energy  region  of  interest,  the  neutron  flight 
times  are  in  the  nanosecond  region,  very  short  ion  bursts  of  high 
intensity  are  needed.  To  satisfy  this  requirement  both  pre-acceleration 
pulsing  and  post-acceleration  compression  have  been  employed  at  the 
University  of  Alberta  5C5  MeV  Van  de  Graaff.  The  deuteron  beam 
delivered  by  the  ion  source  within  the  high  voltage  terminal  is 
chopped  by  means  of  an  RF  of  1MHz  into  pulses  of  10  ns  duration, 
which  after  acceleration  are  compressed  to  0.5  ns  duration  by  means 
of  a  Mobley  magnet  bunching  system.  A  simplified  diagram  of  the  beam 
handling  and  bunching  system  together  with  the  Mobley  magnet  is  shown 
in  fig.  2-1.  The  details  of  the  system  have  been  described  previously 
(Da  66,  Gr  67,  Me  67). 


' 
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Figure  2-1 .  Beam  transport  system. 
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2.2  The  Neutron  Spectrometer 

In  order  to  time  the  neutron  flight  over  an  accurately  known 
distance,  start  and  stop  pulses  are  needed  to  operate  a  Time-to  Ampli¬ 
tude  Converter  (TAC).  The  stop  pulse  originates  from  a  3  cm  long 
cylindrical  capacitor  which  is  30  cm  in  front  of  the  target.  Pulses 
from  the  neutron  detector  start  the  TAC.  The  reason  for  using  the 
detector  pulse  as  the  start  signal  is  that  the  dead  time  of  the  system 
is  decreased.  This  is  because  the  percentage  losses  is  a  product  of 
the  start  rate  and  the  dead  time  (5ysec)  of  the  TAC. 

The  neutron  detector  consist?  of  Ne  218  liquid  scintillator 
in  a  3"  diameter  by  1/2"  long  quartz  container  optically  coupled  to 
an  RCA-4522  photomultiplier.  Using  a  remotely-controlled  cart,  the 
detector  distance  from  the  target  and  its  angle  with  respect  to  the 
incident  beam  can  be  varied  from  1.0  to  6.3  meters  and  from  0  to  150 
degrees  respectively. 

For  normalization  of  angular  distributions  a  fixed-position 
detector  at  2.00  meters  and  30  degrees  serves  as  a  monitor.  It 
consists  of  a  Naton  phosphor  optically  coupled  to  an  RCA-8575  photo¬ 
multiplier  tube.  The  monitor  system  uses  the  same  stop  pulse  as  the 
main  system. 

Fig  2-2  shows  the  block  diagram  of  the  electronics  used  for 
the  time  of  flight  spectrometer.  Signals  from  the  photomultiplier  of 
the  main  detector  are  fed  to  a  Constant  Fraction  Pulse  Height  Trigger 
(Ge  67)  whose  output  is  used  to  start  TAC  #1.  The  signal  from  the 


Figure  2-2. 


A  block  diagram  of  the  electronics  used  for  the 


measurements. 
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beam  pickoff  is  fed  to  a  Fast  Zero-Crossover  Discriminator  (Ge  67a) 
which  stops  TAC  #1  and  #3.  The  output  of  TAC  #1  is  fed  to  an  Analog- 
to-Digi tal  Converter  'A'  (ADCA)  which  feeds  the  SDS-920  on-line  com- 
puter.  TAC  #1  is  strobed  by  the  start  signal  through  the  output  of  a 
Gate  and  Delay  Generator. 

The  pulse  from  dynode  11  (slow  pulse)  is  amplified  and  fed 
to  a  single  channel  analyzer  (SCA)  #1,  which  acts  as  an  integral 
discriminator,  to  set  a  threshold  for  the  neutron  detector.  The  ADCA 
is  gated  by  a  coincidence  between  the  delayed  start  pulse  and  tne 
output  of  SCA  #1.  Since  a  gamma-ray  pulse  is  different  in  shape  from 
that  of  a  neutron  pulse,  an  integration  and  clipping  yields  bipolar 
pulses  whose  zero-crossing  points  are  timed  differently  with  respect 
to  the  leading  edge  of  the  detector  anode  pulse  as  measured  by  TAC 
#2.  Using  SCA  #2  a  window  is  set  on  the  neutron  pulses  which  are 
routed  into  the  second  2048-channel  range  while  the  rest  of  the  pulses 
remain  in  the  first  2048-channels. 

Since  the  conditions  of  high  resolution  and  elimination  of 
y-ray  pulses  can  be  relaxed  in  the  monitor  circuit,  a  simplified 
version  of  the  main  circuit  is  used.  In  order  to  decrease  the  back¬ 
ground  in  the  monitor  Spectrum  a  higher  threshold  is  set  by  SCA  #3. 
This  spectrum  is  stored  in  the  1 1024-channel '  ADC  and  can  be  dumped 
into  the  computer  memory  on  command. 


. 
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CHAPTER  III 

THE  EXPERIMENT  AND  DATA  ANALYSIS 

3. 1  Targets 

Elemental  sulfur  is  not  a  suitable  target  material  since  it 
sublimates  in  vacuum  at  room  temperature.  Several  sulfur  compounds 
such  as  CdS ,  Ag2S  and  Sb2S3  have  been  tried.  It  has  been  found  that 
such  compounds  peel  off  the  backing  before  they  melt,  A  technique 
based  mainly  on  the  method  described  by  D.  D.  Watson  (Wa  66)  has  been 
used  to  prepare  stable  targets. 

Natural  sulfur  (95%  32S,  0.76%  33S,  4.22%  3^S  and  0.02%  36S) 
has  been  evaporated  in  vacuum  onto  a  silver  backing.  The  latter,  with 
the  sulfur  layer  on  it,  was  heated  to  200°C  in  air.  An  immediate  chem¬ 
ical  reaction  occurs  forming  a  very  stable  layer  of  Ag2S.  The  advan¬ 
tage  of  evaporating  sulfur  in  vacuum  rather  than  in  air  is  the  better 
target  uniformity.  Targets  prepared  by  this  method  have  been  proven 
to  be  durable  and  capable  of  withstanding  beam  currents  of  more  than 
1  yamp  at  5.5  MeV  for  many  hours  without  deterioration.  During  the 
experiment  air  cooling  has  been  applied  to  these  targets. 

For  absolute  normalization  of  the  cross  section  a  gas  target 
has  been  used.  At  first  H2S  gas  was  tried  as  a  target  but  it  was 
found  to  dissociate  under  beam  bombardment,  S02  gas  has  been  found 
to  be  stable  and  therefore  was  used  as  a  target. 
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3.2  Calibration  of  the  Machine  Energy 

The  machine  energy  has  been  calibrated  using  Li7(p,n)7Be, 
--CCpjn)1  N  and  i9F(p ,n) i9Ne  thresholds,  The  estimated  error  in  the 
calibration  is  ^  +5  KeV. 

3.3  Calibration  of  the  Time  Spectrum 

The  time  calibration  has  been  done  by  using  the  flight  time 
and  peak  position  of  the  gamma  ray  group  and  finding  the  slope  of  the 
calibration  line  by  the  electronics  setup  of  fig  3-1.  Output  pulses 
from  the  TAC  will  occur  whenever  the  start-stop  time  difference  is  a 
multiple  integer  of  the  period  of  the  oscillator  pulse.  By  measuring 
the  average  spacing  of  the  peaks  (a)  in  channels,  the  slope  of  the 
time  calibration  line  has  been  calculated  from  the  relation: 

Slope  =  10 3/ ( a c f )  ns/channel 
where  f  is  the  oscillator  frequency  (mc/sec). 

3.4  Detector  Efficiency 

The  neutron  detector  efficiency  was  calculated  using  a 
computer  code  developed  by  T.  B.  Grandy  (Gr  67).  The  calculations 
have  been  compared  with  detector  efficiencies  derived  from  experimental 
data  (Bu  67),  and  are  considered  to  be  accurate  to  within  8%  at  the 
neutron  energies  encountered  in  this  experiment. 


. 
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Figure  3-1 


Time-Calibration  Electronics. 
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Since  the  magnitude  of  the  efficiency  is  dependent  upon  the 
lower  threshold  as  determined  by  SCA  #1,  the  cutoff  neutron  energy  must 
be  measured.  The  following  method  has  been  used. 

Using  the  electronics  setup  of  fig  3-2,  a  gamma  ray  spectrum 
from  a  z2Na  source  has  been  obtained.  The  lower  threshold  in  electron 
energy  is  given  by: 

Lower  threshold  =  R.EC 

where  R  is  the  ratio  between  the  cutoff  channel  and  the  channel  corres¬ 
ponding  to  2/3  of  the  height  of  the  0.511  keV  Compton  edge  (fig  3-3). 

Ec  =  340  keV,  is  the  maximum  Compton  electron  energy.  In  order  to  get 
the  lower  threshold  in  terms  of  proton  recoil  energy,  a  response  curve 
of  Ne  218  scintillator  to  electrons  and  protons  is  needed.  Since  such 
a  curve  is  not  avai  1  able ,  the  data  published  by  Smith  ejt  al(Sm  68)  for  Ne  213, 

which  is  quite  similar  to  Ne  218,  have  been  used.  The  overall  error 

in  the  detector  efficiency  has  been  estimated  to  be  10%. 

3.5  The  32S(d,n)33Cl  Reaction 

Using  the  beam  transport  system  and  the  neutron  spectrometer 
(discussed  in  Chapter  II),  the  32S(d,n)33Cl  reaction  has  been  studied 
at  deuteron  energies  of  4.70  and  5.50  MeV,  and  flight  distance  of 
6.325  meters.  An  overall  time  resolution  (FWHM)  of  1.2  nsec  was 
obtained.  With  a  flight  path  of  6.325m,  this  implies  an  energy  resolu¬ 
tion  of  68  keV  for  the  5.5  MeV  neutrons  and  15  keV  for  2  MeV  neutrons. 


i  &  .  •  • 


Figure  3-2.  22Na  Gamma  ray  spectrum. 


Figure  3-3. 


A  block  diagram  of  the  electronics  used  for  measuring 
the  cut-off  neutron  energy. 
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In  order  to  achieve  good  statistics  for  the  weakly  excited 
levels  the  reaction  has  been  repeated  at  a  shorter  distance  of  2.957 
meters  and  4.54  MeV  deuteron  energy.  In  each  experiment  spectra  have 
been  taken  in  the  angular  range  of  0°  to  140°.  Fig  3-4  shows  a 
typical  time  of  flight  spectrum. 

3.6  Excitation  Energies 

The  peak  positions  of  the  various  neutron  groups  were 
obtained  using  a  peak-fitting  program  (Te  66).  For  the  spectra  taken 
at  4.70  and  5.50  MeV  deuteron  energies,  the  program  failed  to  fit  the 
peaks  corresponding  to  the  1.99  and  2.35  MeV  states  due  to  their  small 
peak-to-background  ratios.  Therefore,  their  peak  positions  were 
obtained  by  hand-fitting. 

The  deuteron  energy,  the  time  calibration  line  and  the  peak 
positions  formed  the  input  of  a  least  square  mass-and-Q  program.  This 
program  determines  not  only  the  most  probable  Q-value  of  a  state  but 
the  most  probable  mass  for  the  residual  nucleus.  The  latter  was  very 
useful  in  assigning  levels  at  1.993,  2.351  and  2.688  MeV  to  33C1 
nucleus.  Table  3-1  shows  the  output  of  the  mass-and-Q  program  for  the 
4.54  MeV  data. 

The  average  excitation  energies  of  3  Cl  levels  derived  from 
all  the  data  of  all  the  runs  of  the  three  angular  distributions  are 
listed  in  Table  3-2.  These  excitation  energies  are  based  upon  the 
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Figure  3-4 


A  neutron  time-of-flight  spectrum  for  32S(d,n)33  Cl 
reaction  at  =  4.70  MeV. 
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ground  state  Q-value  of  62  +_  9  keV  obtained  in  the  present  work.  A 
Q-value  of  66  +  12  keV  is  obtained  from  the  mass  difference  (En  67). 
the  energy  loss  of  the  deuteron  beam  in  the  target  has  been  taken  into 
account.  The  errors  were  calculated  from  the  estimated  uncertainties 
in  the  peak  positions  and  time  calibration. 


TABLE  3-1 


Best  Fit  for  the 
Q-value  (keV) 

Best  Fit  for  the 
Residual  Mass  (MU) 

72  +  16 

31.58  +  1.0 

-743  +  11 

31.49  +  0.9 

-1936  +  11 

32.24  +  1.0 

-2304  +  8 

32.25  +  0.9 

-2636  +  4 

31.87  +  0.6 

-2799  +  3 

32.11  +  0.6 

For  comparison,  the  excitation  energies  obtained  previously 
are  also  tabulated  (Table  3-2).  The  results  obtained  in  the  present 
work  agree  fairly  well  with  those  of  Moss  (Mo  67)  and  Prosser  and 
Gordon  (Pr  67)  from  (3He,d)  and  (p,y)  reactions.  However,  the  follow¬ 
ing  discrepencies  appear  between  our  results  and  those  obtained  in  the 
previous  (d,n)  work: 

(i)  We  do  not  observe  either  the  state  at  2.11  MeV 

reported  by  Macefield  and  Towle  (Ma  60)  and  Mubarak- 
mand  and  Macefield  (Mu  67)  or  the  level  at  2.53  MeV 
reported  by  both  of  them  and  by  Middleton  et  al . 
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-A  weak  indication  for  a  peak  was  observed  at  certain  angles. 

-The  excitation  energy  of  this  level  is  obtained  from  36Ar(p,a)33Cl  reaction. 
-The  excitation  energy  is  that  of  Mo  67. 
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(Mi  53).  On  the  other  hand  we  do  observe  states  at 
lc993,  2.351  and  2.688  MeV  which  they  did  not  see 
except  for  a  weak  level  at  1 » 98  MeV  reported  by 
Middleton  et  al_. 

(ii)  In  Grandy  et  aj_.  (Gr  68)  experiment  on  a  42Ca  target 
which  was  contaminated  with  32S,  levels  at  2.12,  2.56, 
2.63  and  2.69  MeV  have  been  identified.  Only  the  2.69 
MeV  state  has  been  confirmed  in  the  present  experiment. 
Comparing  the  peak-to-background  ratio  of  the  spectra 
obtained  in  the  present  work  and  that  of  Grandy  et  al . , 
it  is  found  that  they  are  the  same.  Thus  the  extra 
peaks  observed  by  Grandy  et_  a_l_.  must  have  been  due  to 
some  unknown  contaminant  other  than  32S. 

3.7  Cross  Sections 

Using  plot-and-print  program  (Da  69)  the  background  in  the 
raw  spectra  was  calculated  using  a  parabolic  shape.  Peak  areas  of  the 
various  neutron  groups  were  extracted.  However,  in  the  5.5  MeV  data, 
the  areas  of  the  1.99  MeV  and  2.35  MeV  peaks  could  not  be  obtained 
due  to  the  small  peak-to-background  ratio.  Hand- fitting  has  been 
used  to  obtain  the  areas  of  these  peaks  in  the  4.54  and  4.70  MeV 
data. 

The  neutrons  which  are  misidentified  as  gamma  rays  appear 
as  peaks  in  the  gamma  ray  spectrum.  These  peaks  were  summed  and  added 
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to  the  corresponding  areas  in  the  neutron  part  of  the  spectrum. 

The  contaminant  peaks  due  to  the  presence  of  12C  and  160, 
which  interfere  with  those  of  33C1  at  certain  angles,  were  subtracted. 
The  results  of  N.  Davison  (Da  68)  on  17F  states  and  those  of  A.  W. 
Obst  (Ob  66)  and  C.  E.  Hollandsworth  et  aj_.  (Ho  66)  on  :  3N  states 
have  been  used  for  this  purpose. 

The  peak  areas  obtained  were  normalized  to  the  area  of  the 
ground  state  peak  of  33C1  in  the  monitor  spectrum. 

The  relative  corrected  areas  were  calculated  using  the 
formula: -- 


Relative  Corrected  Area 


tm  _2 

Am  t  e 


where  A  = 
\ 

t  = 


area  of  the  neutron  peak  in  the  neutron  spectrum 

area  of  the  corresponding  neutron  peak  in  the  gamma 
spectrum 

monitor  system  live  time 
main  detector  live  time 

absolute  efficiency  of  the  neutron  detector 


For  the  data  taken  at  4.54  and  4.70  MeV  deuteron  energy  the 
absolute  normalization  of  cross  sections  has  been  carried  out  using  a 
gas  target.  At  first  the  D(d,n)3He  experiment  was  done  to  check  the 
charge  collection  in  the  gas  cell.  A  value  of  67.5  +_  7  mb/sr  was 


i  ; 
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obtained  for  the  cross  section  at  0  and  4-50  MeV  deuteron  energy. 

This  value  agrees  fairly  well  with  the  published  value  of  72  0  +  5 

mb/sr  at  a  deuteron  energy  of  4.60  MeV  (St  56). 

A  run  at  0°  with  S02  gas  at  a  pressure  of  16.45  mm  Hg  was 
used  to  calculate  the  absolute  differential  cross  section.  The  area 
of  the  0-81  MeV  peak  was  extracted  and  corrected  for  dead  time.  The 
corresponding  absolute  differential  cross  section  was  calculated 
using  the  formula: 

(do/dn)  =  1.6  x  1014* (A/QNfte )  mb/sr 

A  =  area  of  the  peak  after  corrections  for  detector  dead 
time  have  been  made 

Q  =  charge  collected  (uC) 

N  =  number  of  32S  nuclei/cm2  present  in  the  gas  cell 

n  =  solid  angle  subtended  by  the  detector  -  1.51  x  10“4 

sr  at  6.3  meters 

e  =  absolute  efficiency  of  the  neutron  detector 

As  a  second  check  on  the  absolute  cross  sections  the  target 
thickness  has  been  measured  Independently  by  determining  the  energy 
width  at  half  maximum  of  the  2.69  MeV  peak.  The  results  of  this 
method  agree  to  within  5%  with  that  of  the  gas  target  method*  For  the 
5,5  MeV  runs,  the  absolute  normalization  has  been  done  using  the  energy 
width  method. 

The  absolute  cross  sections  for  the  states  of  'Cl  are  given 
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in  Tables  3-3  to  3-8, 

The  errors  noted  are  those  due  to  counting  statistics. 

These  errors  are  quite  large  for  the  1.99  MeV  state  due  to  the  small 
peak-to-background  ratio.  The  imperfect  extraction  of  the  2,35  MeV 
peak  from  that  of  Fq  5  was  the  major  source  of  error  in  the  cross 
sections  of  this  state.  Not  shown  are  the  errors  arising  from  the 
detector  efficiency  (10%)  and  the  absolute  normalization  (11%).  Thus 
the  absolute  cross  sections  will  have  an  additional  error  of  15%r 
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DIFFERENTIAL  CROSS  SECTION  FOR  THE  GROUND  STATE 


Ed 

=  4.70  MeV 

Ed  =  5.50  MeV 

c.m.  Angle 

(d0/dn)c.m. 

Error 

c.m.  Angle 

(do/dfi)„ 

'  'c.m. 

Error 

deg* 

mb/sr 

mb/sr 

deg. 

mb/sr 

mb/sr 

0.0 

1.98 

0.07 

0,0 

0.83 

0.06 

5.2 

1,97 

0.10 

5.2 

1.00 

0.07 

10.4 

2.47 

0.08 

10.4 

1.49 

0.06 

15.7 

3,01 

0.13 

20.9 

2.66 

0.11 

20.9 

3.65 

0.10 

31.3 

4.10 

0.14 

31.3 

4.46 

0.13 

41.6 

3.77 

0.12 

41.6 

3.74 

0.10 

51.9 

2.86 

0.12 

51.9 

2.36 

0.08 

62.2 

1.53 

0.06 

62.2 

1.16 

0.05 

72.4 

1.01 

0.05 

72.4 

1.22 

0.06 

82.5 

1.01 

0.07 

82.5 

1,25 

0.05 

92.5 

1.07 

0.04 

102.5 

1.53 

0.05 

112.4 

1,07 

0,07 

122.2 

1 . 40 

0.04 

131.9 

0.82 

0.06 

141.6 

1,20 

0.04 

143.4 

0o  72 

0.09 

TABLE  3-3 
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DIFFERENTIAL  CROSS  SECTION  FOR  THE  0.812  MeV  STATE 


Ed  = 

4.70  MeV 

Ed  =  5,50  MeV 

c,m,  Angle 
deg. 

(da/dn)cm_ 

mb/sr 

Error 

mb/sr 

c,m.  Angle 
deg. 

(d3/d8)c.m. 

mb/sr 

Error 

mb/sr 

0,0 

23.59 

0,52 

0.0 

25.04 

0.90 

5.2 

21.90 

0,49 

5.2 

22.84 

0.82 

10.5 

16.63 

0.37 

10.5 

16.18 

0.58 

15,8 

9,22 

0,26 

20.9 

4.48 

0.16 

21.0 

6,27 

0,17 

31.4 

0.28 

0.02 

31.4 

1.38 

0,07 

41.8 

0.89 

0.07 

41.8 

1.57 

0,07 

52.1 

2.07 

0.08 

52.1 

2. 85 

0.10 

62.4 

2.71 

0,07 

62.4 

2.71 

0,07 

82.7 

1.20 

0.08 

72  c  6 

2. 49 

0.09 

92.7 

0.49 

0.02 

82.7 

1.73 

0.05 

112.6 

0.26 

0.03 

102.7 

0.85 

0.03 

132.1 

0.50 

0.06 

122.4 

0.78 

0.23 

143.5 

0.78 

0.05 

141.8 

1.11 

0.03 

TABLE  3-4 


24 


DIFFERENTIAL  CROSS  SECTION  FOR  THE  1.993  MeV  STATE 


Ed  =  4.54  MeV 

Ed  =  4.70  MeV 

c.m.  Angle 
deg. 

(do/d£5)c.m. 

mb/sr 

Error 

mb/sr 

c.m.  Angle 
deg. 

<d°/dn>c.m. 

mb/sr 

Error 

mb/sr 

0,0 

0.52 

0.05 

0.0 

0a53 

0.08 

5.3 

0.55 

0.06 

5.3 

0.54 

0.07 

10.6 

0.60 

0.05 

10.6 

0.64 

0.08 

15.9 

Oc  63 

0.05 

15.9 

0.39 

0.06 

21.2 

0.64 

0.06 

21.2 

0.56 

0.10 

26.4 

0.62 

0.07 

31.7 

0.57 

0.11 

31.7 

0.52 

0.06 

42.2 

0,58 

0.09 

37.0 

0.48 

0.06 

52,6 

0,62 

0.08 

42.2 

0.40 

0.05 

73.2 

0.46 

0.08 

52.6 

0.40 

0.04 

83.3 

0,41 

0.06 

63.0 

0.35 

0,04 

103.4 

0.39 

0.03 

83.4 

0.51 

0.05 

123.0 

0.38 

0.03 

103.4 

0.34 

0,04 

142.2 

0.38 

0.05 

123.0 

0.47 

0,04 

TABLE  3-5 
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DIFFERENTIAL  CROSS  SECTION  FOR  THE  2.351  MeV  STATE 


Ed  =  4.54  MeV 

Ed  =  4.70  MeV 

c.m.  Angle 

(do/d8)c.m. 

Error 

c.m.  Angle 

(do/dn)c.m. 

Error 

deg. 

mb/sr 

mb/sr 

deg. 

mb/sr 

mb/sr 

0.0 

0.70 

0,10 

0.0 

o,5o; 

0.08 

5.3 

0.72 

0.07 

5,3 

0.59 

0.09 

10.6 

0.97 

0.07 

10.6 

0.73 

0.09 

16.0 

0.97 

0.06 

0*47 

&?e 5 

21.3 

1.06 

0.08 

21.2 

O.J6o.?2 

0.02 

26.6 

1.22 

0.07 

31.8 

0.87 

0.12 

31.9 

1.16 

0.09 

42.3 

0.92 

0.09 

37.1 

0.94 

0.06 

52,8 

0.68 

0.09 

42.4 

0.82 

0,05 

63.1 

0.34 

0.07 

52.8 

0.27 

0,02 

73.4 

0.27 

0.06 

63.2 

0.21 

0.01 

83,6 

0.17  . 

0.04 

83.6 

0.43 

0.06 

103.6 

0,30 

0.04 

103.6 

0.41 

0,05 

123.1 

0.43 

0.04 

123.2 

0.26 

0.03 

142,3 

0.42 

0.06 

142.4 

0.16 

0.04 

TABLE  3-6 
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DIFFERENTIAL  CROSS  SECTION  FOR  THE  2.688  MeV  STATE 


Ed  =  4.70  MeV 

Ed  =  5.50  MeV 

c.m.  Angle 
deg. 

(do/dn)c.m. 

mb/sr 

Error 

mb/sr 

c.m.  Angle 
deg. 

(da/dft) 

c.m. 

mb/sr 

Error 

mb/sr 

0.00 

2.67 

0.07 

0.0 

3.82 

0.17 

5.3 

2.57 

0.07 

5.3 

3.73 

0.13 

10.7 

2.89 

0.08 

10.6 

3.64 

0.16 

16.0 

2.07 

0.11 

21.2 

3.41 

0.15 

21.4 

2.52 

0.11 

31.8 

3.41 

0.15 

32.0 

2.92 

0.11 

42.3 

3.11 

0.10 

42.6 

2.88 

0.10 

52.8 

3.04 

0.13 

53.0 

2.85 

0.10 

63.1 

2.45 

0.08 

63.4 

2.62 

0.07 

73.4 

2.15 

0.09 

73.7 

2.33 

0.08 

83.5 

1.79 

0.11 

83.9 

1.96 

0.04 

93.6 

1.65 

0.03 

103.9 

1.81 

0.04 

113.4 

1.67 

0.07 

123.4 

1.50 

0.05 

132.8 

1.01 

0.04 

142.5 

1.15 

0.03 

144.0 

0.79 

0.08 

TABLE  3-7 
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DIFFERENTIAL  CROSS  SECTION  FOR  THE  2.851  MeV  STATE 


E 

d  =  4.70  MeV 

Ed  =  5,50  MeV 

c.m.  Angle 

(da/do) 

'c.m. 

Error 

c.m.  Angle 

Error 

deg. 

mb/sr 

mb/sr 

deg. 

mb/sr 

mb/sr 

0.0 

33,15 

0.74  1 

0.0 

40.10 

1,44 

5.4 

32.41 

0,72 

5,3 

39.36 

1,41 

10.7 

34.75 

0.77 

10,6 

39.25 

1.41 

16.1 

21,36 

0.47  . 

21.3 

35,88 

1.29 

21,4 

23.02 

0,51 

31.9 

26.04 

0.93 

32,1 

24,56 

0,61 

42.4 

12.66 

0.42 

42.7 

17.39 

0,38 

52.8 

8.09 

0.34 

53.2 

11.70 

0.26 

63.2 

6.48 

0.27 

63.6 

7.25 

0,16 

73.5 

6.18 

0.25 

73.9 

7.72 

0,17 

83,66 

5.56 

0,23 

84,1 

6.90 

0,12 

93,7 

5,62 

0.23 

104.1 

5,91 

0.10 

113.5 

4.03 

0.16 

123.6 

4.50 

0.11 

132,8 

3.31 

0.17 

142.7 

4.10 

0.10 

144.1 

3.35 

0.17 

TABLE  3-8 


. 

. 

. 

- 


28 


CHAPTER  IV 

THEORETICAL  ANALYSIS  OF  ANGULAR  DISTRIBUTIONS 
4.1  Background 

The  two  main  mechanisms  through  which  a  nuclear  reaction 
can  proceed  are  the  compound  nucleus  (CN)  and  the  direct  modes.  Both 
may  simultaneously  contribute  to  the  measured  reaction  cross  section. 
Therefore,  to  analyze  the  direct  part,  the  CN  contribution  should  be 
subtracted  from  the  experimental  cross  section. 

Using  a  statistical  approach  Hauser  and  Feshbach  (Ha  52) 
have  obtained  an  analytic  expression  for  the  differential  cross 
section  of  a  reaction  when  it  proceeds  via  CN  mode.  Their  method  is 
known  as  Hauser-Feshbach  (HF)  theory. 

Direct  reactions  have  been  successfully  explained  by  the 
Distorted  Wave  Born  Approximation  theory  (DWBA).  It  refers  to  the 
method  of  calculating  the  differential  cross  section  of  stripping 
using  the  optical  model  wave  functions.  A  detailed  mathematical 
treatment  of  the  DWBA  theory  has  been  developed  by  a  number  of  authors. 
Excellent  among  these  are  those  given  by  Tobocman  (To  61)  and  G.  R. 
Satchler  (Sa  64) . 

The  greatest  significance  of  the  analysis  of  the  direct 
part  is  that  it  can  yield  the  orbital  angular  momentum  l  of  the  cap¬ 
tured  particle.  The  following  model -independent  selection  rules  can 
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be  used  to  determine  the  change  of  angular  momentum  and  parity  in  a 
stripping  process: 

I  ( I  Ji‘Jf  I  -  Jf  <  Ji  +  l  +  \ 

=  (-1)^ 

where  J.. ,  J^>  and  are  the  initial  and  final  spins  and  parities. 

If  the  spin  of  the  target  nucleus  is  0+,  the  final  state  will  have 

1  7 

spin  =  l  ±  2  and  parity  =  (-1)  . 

The  l- value  of  the  captured  particle  is  obtained  from  fit¬ 
ting  the  direct  part  of  the  experimental  angular  distribution  with  the 
DWBA  predictions.  In  doing  so,  the  ordinates  of  the  theoretical  curves 
have  to  be  adjusted  so  as  to  make  them  match  the  data  at  the  principal 
maximum.  The  DWBA  calculations  are  based  on  the  assumption  that  the 
captured  particle  goes  into  a  single-particle  state,  i.e.,  a  shell 
model  state.  If  this  is  not  the  case,  the  experimental  cross-section, 
of  course,  will  be  lower  than  the  calculated  one.  In  general,  one  can 
wri  te : 


(da/dft)exp  =  S(dc/dft)th 

where  the  factor  S, cal  led  the  spectroscopic  factor,  measures  the  ratio 
of  the  observed  cross  section  to  the  cross  section  expected  for  a 
single-particle  state.  Put  differently,  it  measures  the  percentage 
contribution  of  a  single-particle  state  or  states  with  a  given  angular 
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momentum  l  to  the  total  wave  function  describing  the  residual  state. 


4. 2  Analysis  of  Angular  Distributions 

For  the  reaction  studied  in  this  thesis  a  program  written 
by  W.  R.  Smith  (Sm  65)  and  modified  by  N.  Davison  (Da  68)  to  include 
the  level  density  formulae  has  been  used  to  calculate  the  CN  contri¬ 
bution.  The  resultant  CN  cross  sections  for  the  (d,n)  channel  were 
found  to  be  too  small  (fig.  4-3)  to  affect  the  experimental  cross 
sections  for  any  level  of  33C1.  This  can  be  attributed  to  the  high 
level  density  in  the  (d,p)  channel. 

As  a  result,  the  experimental  angular  distributions  were 
only  compared  to  the  DWBA  predictions.  A  DWBA  code  written  by  P.  D. 
Kunz  (Ku  67)  was  used  to  calculate  the  32S(d,n)33Cl  differential 
cross  sections.  For  the  entrance  and  exit  channels,  optical 
potentials  of  the  type: 

U(r)  =  Uc(r)  -  V  f(r,rQr,ar) 


+  4iai  W^rCftr.r^)] 

♦A'  V.  ■ 


muc 


s  r 


os  s 


-► 

a 


were  used,  where  the  Coulomb  potential  U  is  that  of  a  uniformly 
charged  sphere  of  radius 

l 
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The  form  factors  for  the  other  potentials  are  taken  to  have  the  Saxon- 
Woods  shape 

1 

f(r,rQ,a)  =  {1  +  exp  [ (r  -  rQ  AT)/a]}_1 

V,rQ  and  a  are  the  strength,  radius  and  diffussness  of  the  potential. 
The  subscripts  r,  i  and  s  refer  to  the  real,  imaginary  and  spin-orbit 
potentials,  respectively.  The  wave  function  of  the  captured  proton 
was  obtained  by  solving  Schrodinger  equation  with  the  potential 

U  (r)  =  U  (r)  -  V  {f(r,r  ,a  ) 
p\  /  c  v  or  r 

-  *s0  (fi/2mpc)2  f  [f(r.ros.as)]  t  '  ?> 

where  A$o  determines  the  strength  of  the  spin-orbit  potential.  For  a 
given  binding  energy,  the  computer  code  adjusts  V  until  the  solution 
gives  the  same  binding  energy. 

The  code  DWUCK  containes  the  options  for  finite  range  effect 
and  non-locality  of  potentials.  The  first  is  approximated  by  incor¬ 
porating  in  the  zero  range  approximation  the  radial  correction  factor: 

A(r)  =  {1  +  R2(2mpmn/fi2md)  (Ud(r)  -  Up(r)  -  B^}"1 

where  m  .,  m  and  m  are  the  masses  of  the  deuteron,  proton  and 
d  p  n 

neutron  respectively,  and  Bd  is  the  binding  energy  of  the  deuteron. 

R  is  the  finite  range  parameter.  The  non-locality  of  the 
potentials  is  taken  into  account  by  multiplying  the  wave 
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functions  generated  by  the  equivalent  local  potentials  by  factors  of 
the  form: 

i 

F(r)  =  {1  -  B2(m/2fi2)U(r)}"r 

where  e  is  the  non-locality  parameter. 

The  initial  deuteron  optical  model  parameter  used  were 
those  obtained  by  Perey  and  Perey  (Pe  66)  from  the  elastic  scatter¬ 
ing  of  11.8  MeV  deuterons  on  natural  sulfur.  Since  these  parameters 
are  obtained  at  considerably  higher  energy  than  that  at  which  the 
present  work  was  carried  out,  the  depth  of  the  imaginary  well  has 
been  decreased  to  reduce  the  probability  of  absorption  of  the  deuteron 
by  the  target  nucleus. 

Neutron  parameters  for  the  outgoing  channel  are  those  of 
Perey  and  Buck  as  compiled  by  Rosen  (Ro  66). 

Table  4-1  gives  the  optical  model  parameters  used  in  both 
HF  and  DWBA  calculations. 

Figures  4-1  to  4-6  show  the  experimental  data  and  the  DWBA 
predictions  for  the  differential  cross  section  for  33C1  levels.  The 
DWBA  predictions  were  calculated  with  local  potentials  (e  =  0),  a 
lower  cutoff  of  zero  fm  in  the  radial  integral  and  a  value  of  R  = 

0.621  fm  for  finite  range  correction. 
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OPTICAL  MODEL  PARAMETERS  USED  IN  DWBA  AND  HF  CALCULATIONS 


Parameter 

fa) 

Deuteronv  ' 

Deuteron^ 

Captured 

Proton 

Neutron 

Proton^ 

Alpha(c) 

V(MeV) 

154.6 

154.6 

Adjusted  by 
the  computer 
code 

47.0 

54.1 

200,0 

ror(fm) 

0.785 

o 

o 

-'J 

00 

cn 

1.25 

1.27 

1,25 

1.69 

ar(fm) 

0.953 

0.953 

0.65 

0.86 

0,65 

0,576 

W(MeV) 

14.6 

9.0 

9.6 

13.5 

9,0 

~s 

o 

? 

1.567 

1.567 

1.27 

1.25 

1,69 

ai  (fm) 

0.534 

0.534 

0,47 

0.47 

0,576 

Vs(MeV) 

5.91 

5.91 

7.2 

7.5 

0,0 

r  (fm) 
os 

0.785 

1.25 

1.27 

as(fm) 

0,953 

0o65 

0.66 

roc(fm) 

1.3 

1.25 

C n 

O 

28 

1 

(a)  Optical  Model  Parameters  of  Perey  and  Perey  (Pe  66). 

(b)  Optical  Model  Parameters  used  in  the  present  work. 

(c)  The  proton  and  alpha  particle  optical  model  parameters 
were  used  only  In  the  Hauser-Feshbach  Calculation. 
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Figure  4- 


Figure  4-2 


Figure  4-3 


Figure  4-4 


Figure  4-5 


Figure  4-6 


1.  Angular  distribution  of  the  neutrons  leading  the  ground 
state  of  33C1  for  Ed=4.70  and  5.50  MeV.  The  error  bars 
indicate  the  statistical  errors.  The  solid  curves  are 
the  DWBA  predictions. 

.  Angular  distribution  of  the  neutrons  leading  to  the 
0.812  MeV  state  of  33C1  for  Ed=4.70  and  5.50  MeV  (see 
also  caption  under  figure  4-1. 

Angular  distributions  of  the  neutrons  leading  to  the 
1.993  MeV  state  of  33C1  for  Ed=4.54  and  4.70  MeV.  The 
dashed  curve  is  the  predictions  of  the  HF  calculations 
for  Ed=4.54  assuming  JTT=5/2+. 

.  Angular  distribution  of  the  neutrons  leading  to  the 

2.351  MeV  state  of  33C1  for  Ed=4.54  and  4.70  MeV  (see  also 
caption  under  figure  4-1). 

Angular  distribution  of  the  neutrons  leading  to  the 
2.688  MeV  state  of  33C1  for  Ed=4.70  and  5.50  MeV  (see 
also  caption  under  figure  4-1). 

Angular  distribution  of  neutrons  leading  to  the  2.851 
MeV  state  of  33C1  for  Ed=4.70  and  5.50  MeV.  Points 
without  error  bars  indicate  that  the  errors  are  -  the 
size  of  the  points  (see  also  caption  under  figure  4-1). 
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4.3  Discussion 

Satisfactory  fits  have  been  obtained  assuming  an  l  =  2  and 
Z  =  0  for  the  ground  and  the  0.81  MeV  respectively.  It  is  therefore 
possible  to  confirm  that  the  spin-parity  (J77)  is  3/2+  and  l/2+ 
respectively,  in  agreement  with  previous  work. 

The  1.99  MeV  state  is  weakly  populated  in  this  reaction. 

The  angular  distributions  obtained  show  no  stripping  pattern  so 
that  we  are  unable  to  resolve  the  ambiguity  of  J  =  5/2  or  7/2  for 
this  level  (Pr  67) . 

All  the  levels  above  the  1.99  MeV  state  are  unbound  to 
proton  emission.  However  higher  levels  of  33C1  can  be  considered 
as  "partially  bound"  owing  to  the  presence  of  the  combined  Coulomb 
and  centrefugal  barrier.  The  code,  DWUCK,  used  in  this  work  includes 
the  option  of  analyzing  the  unbound  "bound"  states. 

The  2.35  MeV  level  is  also  weakly  populated  in  this  reaction 
and  is  expected  to  have  an  Z  -  2  in  analogy  to  the  2.313  MeV  state  in 
33S.  An  Z  =  1  and  Z  =  2  have  been  tried.  The  presence  of  17F  first 
excited  state  which  has  a  strong  Z  =  0  transition  made  a  definite 
assignment  of  this  level  impossible. 

In  the  (3He,d)  work  of  A.  Graue  (Gr  66)  Z-values  of  2  and 
3  were  tried  for  the  level  at  2.69  MeV.  No  definite  Z- value  could 
be  obtained  since  the  angular  distributions  predicted  by  the  DWBA 
theory  used  for  the  (3He,d)  work  are  not  much  different.  In  the 
present  work  we  find  that  the  predicted  angular  distributions  for 
l  =  2  and  3  are  quite  different  for  the  (d,n)  reaction.  It  is 
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clear  from  fig.  4-5  that  the  2.69  MeV  level  shows  an  angular  distribution 
consistent  with  7  =  3  and  a  level  spin  of  5/2"  or  7/2'.  The  2.69  MeV 
state  is  very  likely  to  be  the  1 f^2  state  expected  in  this  range  of 
excitation  energy.  The  theoretical  work  of  Watson  and  Lee  (Wa  67) 
predicts  a  lfy  2  state  in  33C1  at  2.75  MeV  excitation. 

The  2.85  MeV  state  is  known  from  gamma  ray  work  (V a  58,  Pr  67) 
to  be  a  doublet  with  spins  5/2+  and  3/2".  The  angular  distributions 
obtained  in  the  present  work  shows  that  the  7  =  1  member  is  much  stronger 
than  that  of  the  7=2.  This  result  agrees  with  that  obtained  by  Graue. 

The  absolute  spectroscopic  factors  extracted  using  the  data 
at  4.70  and  5.50  MeV  are  given  in  table  4-2.  The  values  of  S  obtained 
at  the  two  energies  agree  within  the  indicated  uncertainties. 

The  experimentally  determined  spectroscopic  factors  depend  on 
the  cross  section  predicted  by  the  DWBA  theory.  It  is  important  to  find 
out  how  sensitive  the  calculated  cross  section,  and  hence  the  spectrosco¬ 
pic  factors,  are  to  changes  in  some  of  the  optical  model  parameters,  the 
incorporation  of  nonlocality  approximation,  and  the  omission  of  finite 
range  correction.  This  has  been  investigated  and  the  results  are  pre¬ 
sented  in  table  4-3.  S  seems  to  be  insensitive  to  spin-orbit  effects, 
finite  range  correction  and  to  large  changes  in  the  imaginary  well  depth 
for  the  deuteron.  However,  the  addition  of  nonlocality  creates  a  change 
of  <_  27%  in  S  but  it  should  be  noted  that  its  incorporation  is  only  an 
approximation.  We  can  then  see  that  although  the  deuteron  optical  para¬ 
meters  are  obtained  from  data  at  bombarding  energy  much  higher  than 
that  used  in  the  present  experiment,  the  resultant  spectroscopic  factors 
should  not  be  greatly  affected  by  this. 
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Change  in  Optical  Model 

Percentage  change  in  $(£^=4.7  MeV) 

Paramete  rs 

GND 

0.81 

2.69 

2.85 

7=2 

7=0 

7=3 

7=1 

With  WD  =  14.6  MeV 

+6 

+  16 

+  10 

+9 

(62%  change) 

•j* 

With  nonlocal  potentials 

-25 

-27 

-20 

-16 

Without  spin-orbit  for  the 

+5 

+6 

+3 

+4 

deute ron 

Without  spin-orbit  for  the 

+20 

0 

+  16 

+3 

Captured  Proton 

Without  spin-orbit  for  the 

-3 

+2 

+1 

+  1 

neutron 

With  zero  range 

+3 

+3 

+2 

+5 

TABLE  4-3 


Sensitivity  of  the  Spectroscopic  Factors  to  the  Form 
and  Parameters  of  the  DWBA  Calculation 


t  en  =  3p  =  0.85  fm. 


3d  =  On 54  fm. 
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For  comparison,  the  values  of  S  obtained  by  Graue  (Gr  66) 
using  the  (3He,d)  reaction  are  also  shown  in  table  4-2.  Although 
these  values  are  somewhat  different  from  what  we  have  measured  the 
general  trend  of  both  sets  is  in  agreement .  Since  the  spin  of  the 
2.69  MeV  state  is  not  known  Graue  reported  the  value  of  the  strength 
function  ( 2 J+l ) S  rather  than  the  spectroscopic  factor.  Thus  to  compare 
his  result  for  this  level  with  that  of  the  present  work  we  have 
calculated  S  assuming  =  7/2"  , 


4.4  Shell  Model  Description  for  33C1  Nucleus 

Glaudemans  et  a_L  (G1  64)  have  carried  out  shell  model  calcula 
tions  for  levels  of  even  parity  in  nuclei  in  the  mass  range  29  to  40. 

An  inert  28Si  Core  is  assumed  with  a  two-particle  interaction  of  the 
outer  nucleons  in  the  2Sy2  ar|d  ^3/2  s*ie^s°  Mixing  of  all  possible 
configurations  in  these  shells  is  taken  into  account.  The  level  scheme 
predicted  by  these  calculations  for  A- 33  is  shown  in  fig.  4-7.  By  using 
a  delta  function  as  the  effective  two-nucleon  interaction  acting 
only  in  the  nuclear  surface  as  has  been  suggested  by  Green  and 
Moszkowski  (Gr  65)  ,  Glaudemans  et  al_.  (G1  66)  recalculated  the 
excitation  energies,  and  the  new  energy  values  agree  very  well  with 
the  experimental  energies.  However,  the  odd  parity  levels  are  not 
predicted  by  these  calculations.  Erne  (Er  66)  has  performed 
shell  model  calculations  of  levels  of  odd  parity  in  the  region  of  the 
lcL/?  shell  nuclei.  An  inert  i2S  Core  is  assumed  with  a  residual  two- 
particle  interaction  of  the  outer  nucleons.  Only  one  nucleon  is 
considered  to  be  in  the  lf^2  shell,  while  the  others  are  in  the  ld^ 
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Figure  4-7.  Comparison  of  the  theoretical  energy  levels  and 
spins  with  the  experimental  data  on  the  low-lying 
states  in  the  ?3C1  nucleus „ 
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shell.  According  to  these  calculations,  the  If ?/2  state  in  33C1  is 
predicted  at  2 . 9 4+0 . 5  MeV  excitation.  The  level  observed  in  the  present 
work  at  2.69  MeV  excitation  is  most  probably  the  If ^  ^  state. 

It  is  interesting  to  find  out  how  the  shell  model  can  be  used 
to  make  a  quantitative  analysis  of  spectroscopic  factors.  S  is  defined 
theoretically  as  an  overlap  integral  given  by: 

Sh(lj)  =  (A+l)*1  Hy  d(A+l) 

where  $  is  the  wave  function  constructed  by  the  vector  coupling  to 
spin  Jf  of  a  particle  in  the  state  l j  to  the  target  wave  function  of 
spin  .  The  symbol  d(A+l)  denotes  integration  over  the  (A+l)  nucleons 
in  the  final  state.  Thus  by  overlapping  the  wave  functions  predicted 
by  the  shell  model  for  the  initial  and  final  states,  predictions  of 
spectroscopic  factors  can  be  obtained. 

In  the  early  calculations  of  Glaudemans  et_  aj_.  (G 1  64),  the 
wave  functions  of  the  ground  state  of  32S  and  the  predicted  states  of 
33C1  have  the  following  major  configurations:  - 
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1110 

0131 

1031 
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where  the  notation  of  Glaudemans  et  al .  has  been  used.  The  number  in 
brackets  under  the  first  column  refer  to  the  calculated  excitation 
energy  of  the  level. 

Since  the  wave  functions  have  been  determined,  the  spect¬ 
roscopic  factors  for  the  stripping  reaction  32S(d,n)33Cl  could  be 
calculated.  The  general  expressions  for  S  have  been  discussed  by 
Macfarlane  and  French  (Ma  60a).  The  values  of  S  calculated  by 
Glaudemans  et  aj_.  (G1  64)  are  presented  in  table  4-2. 

The  value  of  S  for  the  ground  state  agrees  with  our  result. 
The  predicted  value  of  S  for  the  first  excited  state  is  a  factor  of 
2  higher  than  the  measured  value.  With  the  assumption  that  only 
2s-|^2  anc*  ^3/2  nuc^eons  take  part  in  the  stripping  process,  the 
transition  to  the  second  excited  state  is  impossible  due  to  angular 
momentum  relation.  A  small  value  of  S  has  been  predicted  for  the 
third  excited  state.  The  calculated  spectroscopic  factors  for  the 
second  and  third  excited  states  are  consistent  with  the  observed 
weak  transition  to  the  1.993  and  2.351  MeV  states. 
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CHAPTER  V 
CONCLUSION 

Three  levels  in  33C1  at  1.993,  2.351  and  2.688  MeV  excitation 
has  been  seen  for  the  first  time  in  the  32S(d,n)33Cl  reaction.  These 
results  confirm  those  obtained  recently  using  (3He,d)  and  (p,y) 
reactions  on  32S.  We  do  not  observe  states  at  2.1  and  2.5  MeV  as 
reported  in  previous  (d,n)  work. 

Angular  distribution  study  confirms  the  previous  spin-parity 
assignment  of  3/2+  and  l/2+  for  the  ground  and  0.812  MeV  states, 
respectively.  Also,  it  confirms  that  one  member  of  the  known  doublet 
at  2.85  MeV  observed  in  (p,y)  experiments  has  a  spin-parity  of  3/2”. 

For  the  weakly  populated  level  at  1.993  MeV  no  Z-value  could 
be  obtained  due  to  the  absence  of  stripping  structure  in  its  angular 
distribution.  An  ambiguity  also  remains  of  l- 1  or  2  for  the  level 
at  2.351  MeV  excitation. 

An  1- 3  has  been  obtained  for  the  level  at  2.688  MeV  and 
hence  its  spin  is  limited  to  be  5/2“  or  7/2”. 

The  spectroscopic  factor  predicted  for  the  ground  state 
using  shell  model  calculations  agrees  with  what  we  have  measured. 

A  qualitative  agreement  has  been  obtained  for  the  0.812,  1.993  and 
2.351  MeV  states.  This  suggests  that  33C1  nucleus  can  be  well 
described  by  a  28Si  core  with  outer  nucleons  in  the  2s^2  ancl  ^3/2 


shells . 
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